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' Abstract 

o ' 

■ We ask whether the universe can be a patchwork consisting of 

, distinct regions of matter and antimatter. We demonstrate that, af- 
ter recombination, it is impossible to avoid annihilation near regional 

^ boundaries. We study the dynamics of this process to estimate two of 

Q_|. its signatures: a contribution to the cosmic diffuse 7-ray background 

^ I and a distortion of the cosmic microwave background. The former 

' signal exceeds observational limits unless the matter domain we in- 

c/3 . habit is virtually the entire visible universe. On general grounds, we 

^ ' conclude that a matter-antimatter symmetric universe is empirically 

excluded. 

- ■ 1 Introduction and Outlook 



The laws of physics treat matter and antimatter almost symmetrically, and 
yet the stars, dust and gas in our celestial neighborhood consist exclusively 
of matter. The absence of annihilation radiation from the Virgo cluster 
shows that little antimatter is to be found within ~20 Mpc, the typical size 
of galactic clusters. Furthermore, its absence from X-ray-emitting clusters 
implies that these structures do not contain significant admixtures of matter 
and antimatter. 

*cohen@bu. edu, derujulaOnxth.21 . cern. ch, glashowOphysics . harvard . edu 
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Many cosmologists assume that the local dominance of matter over anti- 
matter persists throughout the entire visible universe. A vast literature at- 
tempts to compute the baryonic asymmetry from first principles. However, 
observational evidence for a universal baryon asymmetry is weak. In this re- 
gard, searches for antimatter in cosmic radiation have been proposed 
Early next century, the AntiMatter Spectrometer (AMS), deployed aboard 
the International Space Station Alpha 0, |^, will search for antimatter in 
space. Its reach is claimed to exceed 150 Mpc The detection of cosmic 
anti-alpha particles would indicate the existence of primordial antimatter; 
the detection of anti-nuclei with Z > 2 would imply the existence of extra- 
galactic anti-stars. 

The possible existence of distant deposits of cosmic antimatter has been 
studied before 0-0. Steigman concluded that observations exclude sig- 
nificant matter-antimatter admixtures in objects ranging in size from planets 
to galactic clusters. Stecker et al. interpreted an alleged shoulder in the 
cosmic diffuse gamma (CDG) spectrum near 1 MeVQ as relic 7-rays from an- 
timatter annihilation. Recently, Dudarewicz and Wolfendale used similar 
arguments to reach a contrary conclusion: that the observed CDG spectrum 
rules out any large antimatter domains. These conflicting results are not 
based on specific dynamics in a consistent cosmology. Our analysis uses cur- 
rent data and avoids ad hoc assumptions concerning a matter-antimatter 
universe. 

We explore the possibility of universal (but not local) matter-antimatter 
symmetry. In what we term the 5 = universe, space is divided into re- 
gions populated exclusively by matter or antimatter. Our conclusions do 
not depend on how this structure evolved, but it is reassuring to have an 
explicit model in mind: consider an inflationary cosmology in which baryon 
(or antibaryon) excesses develop in the manner suggested by Sakharov Q. 
In models with spontaneous CP violation, the Lagrangian may be chosen 
judiciously so that the 'sign' of CP violation (determining whether a local 
baryon or antibaryon excess develops) is randomly and abruptly assigned 
to regions as they emerge from their horizons during inflation. Soon after 
baryogenesis, the domain walls separating matter and antimatter evaporate. 
As regions of matter or antimatter later re-enter their horizons, the i? = 

"'^We refer to cosmic diffuse photons by conventional names according to the current 
photon energy: the night sky, the CDG and the CBR refer to visible, ~ 1 MeV and 
microwave photons, respectively. 
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universe becomes a two-phase distribution. 

Let today's domains be characterized by a size do, such that I/c/q is 
their mean surface-to-volume ratio. Because the existence of anti-galaxies 
within a matter-dominated domain is empirically excluded, we must (and 
can) arrange the distribution of domains to be sharply cut off at sizes smaller 



than (Iq. Explicit inflationary models satisfying these constraints exist |T0 
but are described no further because we find all such models to conflict with 
observation. 

The current domain size do is the only parameter of the B = uni- 
verse crucial to the confrontation of theory with observation. To agree with 
constraints from X-ray-emitting clusters, do must exceed a minimal value, 
~ 20 Mpc. For do = 20 Mpc, the visible universe would consist of ~ 10^ 
domains. We derive a stronger lower limit on do comparable to the current 
size of the visible universe, thereby excluding the B = universe. 

An explicit cosmological model is necessary to estimate the observable 
signals produced by annihilation. We assume a Robertson-Walker universe 
and use fiducial values for the relevant cosmological parameters: critical mass 
density = 1; vanishing cosmological constant fl\ = 0; Hubble constant 
Ho = 75 km/s-Mpc oi h = 0.75; and an average baryon (or antibaryon) 
number density ub = f]n^ with rj = 2 x 10^^°. 

In Section ^ we show that our conclusions are unaffected by other choices 
for Q, Q\ and Hq within their empirically allowed domains. Consequently we 
do not express our results explicitly in terms of these cosmological parame- 
ters. The annihilation signals we study depend linearly on rj. To compute 
lower limits to the signals, we chose rj at the low end of the domain allowed 
by analyses of primordial element abundances. (For a recent discussion, see 
0.) 

In Section ^ we explain why particle-antiparticle annihilation is unavoid- 
able from the time of recombination to the onset of structure formation. 
Following a conservative approach, we consider only those annihilations oc- 
curring during this period. Our analysis involves known principles of particle, 
atomic and plasma physics, but the dynamics of the annihilating fluids (dis- 
cussed in Section ^ is complicated and the considerations required to reach 
our results are elaborate. 

The immediate products of nuclear annihilation are primarily pions (vr"*", 
TT^ and TT^) with similar multiplicities and energy spectra. The end products 
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are energetic photons from vr*^ decay, energetic electrons^ (e"*" and e~) from 
the decay chain vr /i e, and neutrinos. Although they are produced 
at cosmological distances, the annihilation photons and electrons can each 
produce potentially observable signals: 

• The energy carried off by annihilation electrons (about 320 MeV per 
annihilation) affects the CBR spectrum directly (via Compton scatter- 
ing) and indirectly (by heating the medium). The consequent distortion 
of the CBR (discussed in Section H) cannot exceed observational limits. 

• Most of the annihilation photons, although redshifted, are still present 
in the universe. Their flux (computed in Section ^ cannot exceed the 
observed CDG flux. 

Because these signals increase inversely with the domain size, our analysis 
yields a lower limit for do. In fact, we obtain no new constraint from compar- 
ing the expected distortion of the CBR with its measured limits. However, 
the CDG flux produced by annihilation far exceeds the observed flux unless 
(Iq is comparable in size to the visible universe. Thus, the B = universe is 
excluded. 

2 The Era of Unavoidable Annihilation 

What if the matter and antimatter domains are and have always been spa- 
tially separated? If large empty voids lay between them, there would be no 
observable annihilation signals. We now show how the observed uniformity 
of the CBR rules voids out. 

Two events took place at roughly the same time in cosmic history: the 
transition from charged plasma to neutral atoms (recombination) and the 
decoupling of radiation and ordinary matter (last scattering). For our fiducial 
cosmological parameters, these events occurred at a temperature ~0.25 eV 
and at a redshift ur ~ 1100 (we use y = 1 + z = 1/R{t) as a redshift 
parameter, rather than the conventional z). The transition to transparency 
was not instantaneous, but evolved during an interval ± 100 whose half- 
width is ~ 15 Mpc in comoving (current) distance units. Thus, features at 
recombination of comoving size smaller than 15 Mpc cannot be discerned in 
the CBR. 

^Relativistic electrons and positrons behave similarly, and we refer to both as electrons. 
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Large-scale non-uniformities of the matter density, whether dark or bary- 
generate variations of the CBR temperature. Its observed unifor- 



mity (to parts in 10~^) imphes a very uniform density of ordinary matter 
ai y = yji ^ 1100, to within the resolution discussed above. It follows that 
voids between matter and antimatter domains must be smaller than 15 Mpc. 

The baryon density depletion in voids is damped as photons diffuse toward 
less dense regions, dragging matter with them. By recombination, inhomo- 
geneities with current size < 16 Mpc would be destroyed^ by this mechanism 



||13|| . This upper bound coincides with the smallest resolvable structure in the 
CBR. Thus, voids large enough to survive until recombination would have 
been detected. While matter and antimatter regions may have been sepa- 
rated prior to recombination, they must be in immediate contact afterward. 
Thus, in determining the minimal signal of a -B = universe, we do not 
consider annihilations occurring a.t y > yn. 

The mechanism by which the nearly uniform universe at large y evolved 
today's large-scale structures is not well understood. We cannot confidently 
assert what effects this will have on annihilation in a i? = universe. It could 
well be that the collapse of baryonic matter into galaxies and stars quenches 
annihilation unless the collapsing system overlaps a domain boundary, a sit- 
uation we consider shortly. Our conservative estimate of the annihilation 
signal includes matter-antimatter annihilation taking place prior to the red- 
shift ys at which the earliest density fluctuations become large {Sp/p ~ 1). 
We take ys — 20, which is estimated to be the epoch of galactic conden- 



sation and earliest star formation |]I5[. We compute the signals due to 
annihilations taking place during the interval 1100 > y > 20. 

The large-scale density contrast of the visible universe need not coincide 
with the pattern of matter and antimatter domains. A density fluctuation 
beginning to collapse could overlap a domain boundary. Successful collapse 
would yield a structure with a significant mixture of matter and antimatter. 
In this case annihilation would proceed even more rapidly at the onset of 
structure formation. Yet we cannot be confident that such mixed structures 
form. 

In the linear regime {Sp/p <^ 1), the mean annihilation rate is not af- 
fected by density fluctuations. But, what happens as the fluctuations grow? 
If an over-density is to overcome expansion and become a self-gravitating 



■^This assumes that such inhomogeneities are not strictly isothermal, a situation con- 
sidered in Section H. 
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system, it must satisfy the Jeans condition: the sound travel time across the 
object l/vs must be greater than the characteristic free-fall time 1/ \/G p, or 
GP p > Vg. Suppose that equality is approached by an over-density con- 
taining both matter and antimatter. Further contraction increases the an- 
nihilation rate, thus reducing p and driving the system away from collapse. 
Thus, our conservative estimate of the annihilation signal assumes that den- 
sity fluctuations straddling domain boundaries either fail to collapse or form 
separate unmixed structures. 

3 The Matter— Antimatter Encounter 

The B = universe consists of matter and antimatter domains with almost 
identical mass densities that, as we have shown, touch one another from re- 
combination to the onset of structure formation. As annihilation proceeds 
near an interface, a flow develops as new fluid replenishes what is annihi- 
lated. This flow must be analyzed to determine the annihilation rate on 
which our putative signals depend. The analysis involves established and 
well- understood principles of physics, but is complicated by the energy re- 
leased by nuclear annihilation. (We neglect e~^e~ annihilation, whose energy 
release is much smaller.) The processes by which annihilation electrons lose 
energy produce crucial effects on the ambient fluid, as well as a potentially 
observable distortion of the CBR. (High-energy photons from 7r° decay, al- 
though responsible for the CDG signal, have little effect on the medium 
through which they pass.) 

The primary energy-loss mechanism of the annihilation electrons is Comp- 
ton scattering off CBR photons (see Appendix 0). This process up-scatters 
target photons to higher energies. The resultant flux of UV photons heats 
and ionizes ambient matter throughout much of the universe and for all of 
the relevant period. Moreover, the annihilation electrons lose a small portion 
of their initial energies by scattering off ambient electrons in the fluid. This 
process heats the fluid within the electron range, thereby accelerating the 
flow and leading to even more annihilation — a feedback mechanism making 
the matter-antimatter encounter potentially explosive. 

Several length scales characterize the fluid dynamics about a matter- 
antimatter interface. They are: 

• A, the width of the annihilation zone, wherein both matter and anti- 
matter are present; 
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D, the width of the depletion zone, wherein fluid flow toward the an- 
nihilation zone reduces the density; 

L, the width of the reheated zone, wherein electrons produced by an- 
nihilations directly deposit energy into the fluid. This is simply the 
electron range. 



These length scales, computed in Appendix |A.1| and later in this section, are 



shown in Fig. |T| along with a comoving domain size of 20 Mpc and the horizon 
scale. Annihilation takes place in the vicinity of the domain boundary and 
well within the depletion zone, which itself is much shorter than the electron 
range. That is, in the relevant redshift domain: A <^ D <^ L. This distance 
hierarchy lets us treat the flow as one-dimensional. 

Annihilation has a negligible effect on the CBR temperature T^{y), which 
remains as it is in a conventional universe. However, the annihilation debris 



produce and maintain virtually total ionization, as shown in Appendix |A.3 . 
Therefore the annihilating fluid consists of photons, protons, antiprotons, 
electrons and positrons^. The proton and electron number densities coincide, 
except in the narrow annihilation zone. Consequently our analysis may be 
put in terms of the total matter mass density p = m^nf. + rripUp, the total 
fluid momentum density pv, the total fluid pressure p and the total fluid 
energy density e = ethermai + pv^/2. The internal energy density and pressure 
are related as for a non-relativistic ideal gas: ethermai = Sp/2. 

The equations describing the flow of matter are conservation laws for 
particle number (mass in the non-relativistic limit), momentum, and energy. 
They must take account of the following phenomena: 

• The depletion of fluid mass, momentum and energy by nuclear annihi- 
lation. 

• The effect of the CBR on the fluid momentum and energy. 

• The effect of the annihilation products on the fluid momentum and 
energy. 

• The expansion of the universe. 

^We neglect the helium contamination (^^ 7 % by number) and those of larger primordial 
nuclei. 



7 



10 



10" 



10" 



20 Mpc domain 

Horizon 

Electron Range 

Depletion Zone 

Annihilation Zone 



10 



100 



1000 



y=l+z 



Figure 1: Read from the top at large y: the horizon, the look-back size of 
a 20 Mpc domain, the widths of the reheated zone, the depletion zone, and 
the annihilation zone. 



The expansion of the universe is taken into account by expressing the con- 
servation laws in a Roberston- Walker universe [0. The metric is ds^ = 



dt^ - R^{t)dx^, with X a comoving spatial coordinate normal to a domain 
boundary. The remaining effects are dealt with by including appropriate 
source terms in the fluid equations: 

dp 1 dipv) ^ 

dipv) , R 19,9-, ^ A 
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Here Fann = {cann'v)np IS the matter annihilation rate, ax is the Thompson 
cross section and is the CBR energy density. The terms involving 0"^ 
describe the transfer of energy and momentum between the fluid and the CBR 
resulting from Compton scattering. H^, given by Eq. ( PI) in Appendix ^ is 
the rate of change of the energy density of the fluid due to its interactions with 
the annihilation debris. It receives a direct contribution from the annihilation 
electrons, and an indirect one from UV photons up-scattered by Compton 
collisions of these electrons with the CBR. We find that the contribution of 
the electrons dominates within the electron range. Beyond this range, only 
the UV photons contribute to H^. In Eq. (H), we have neglected the small 
contribution by the annihilation debris to the fluid momentum. 

The signals of a 5 = universe — the CDG and a distortion of the CBR — 
are functions of J, the number of annihilations taking place per unit time 
and area orthogonal to the surface of an annihilation zone: 

J = J (0"annfpp) UpUp R dx , (4) 

where the integral extends over a single annihilation zone with % = at its 
mid-point, and Upix) = np{—x)- The width of the annihilation zone A may 
be estimated as A ~ J/ {<Jppv)nl^, where [|1^ appV ~ 6.5 x 10"^"^ cm^s""*^ c/v 



and rioo is the proton density far from the annihilation zone. We must solve 
Eqs. (11^-0) to determine J. 



3.1 A Qualitative Solution 

Because our fluid equations do not admit analytic solutions we begin with 
a qualitative discussion. The value of Fann is always much greater than the 
expansion rate, so that the solutions to Eqs. rapidly reach equilibrium 

in the annihilation zone. Consequently, taking the limit cTann — ^ oo yields a 
good approximation. In this limit, the width of the annihilation zone shrinks 
to zero and the annihilation terms in the fluid equations may be replaced by 
a boundary condition at the domain interface. The rate of annihilation per 
unit surface area is then given by the proton flux at the interface: 



J ~ v{x,t) 



up 



(5) 
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Two effects result from the couplings of the fluid to the CBR. The term 
in Eq. (Q) proportional to (Tt tends to damp the fluid motion. The corre- 
sponding term in Eq. (|^) tends to keep the fluid temperature near Ty. For 
y > 400 these terms dominate, so that the two temperatures are locked to- 
gether, T c:=L T^. The CBR drag on the fluid leads to diffusive motion, and 
we may define a time- dependent diffusion constant: 

The solution to the resulting diffusion-like equations gives an estimate of the 
annihilation rate J: 

J^rUt)^ (7) 
with rioo the proton number density far from the interface. The width of the 



depletion zone is comparable to the diffusion length D ~ yDe-y t. 

For redshifts y < 200 the effects of the CBR on the fluid motion are 
negligible and we may ignore terms proportional to ctt- In this case, which 
we refer to as 'hydrodynamic', the motion is controlled by pressure gradients 
and the fluid flows at a substantial fraction of the speed of sound. The 
resulting equations are those describing a gas expanding into a semi-infinite 
vacuum in the presence of an energy source. An analytic solution exists for 
Hf = 0. In this case the annihilation rate J is 



J = C riooit) foo(t) = C riooit) 



with Voo the speed of sound and Too the fluid temperature T = p/ irip + rie) far 
from the annihilation zone. The coefficient of proportionality is[| C = (3/4)^. 
The width of the depletion zone in this case is comparable to the sound-travel 
distance D ~ _R(t) dt' v^o/ R- 

In the intermediate region, 200 < y < 400, neither of the above approxi- 
mations give a quantitatively accurate picture of the fluid motion. 



^This is the adiabatic solution. For 100 <y< 200 the process is more nearly isothermal. 
The corresponding value of C is 1/e. 
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3.2 The Numerical Solution 



We have integrated Eqs. (|l|-§|) numerically to determine the fluid temper- 
ature T and the annihilation rate per unit surface area J near a domain 
boundary. The diffusive nature of the solution at large y has a welcome 
consequence: all memory of the initial conditions is lost as the fluid evolves. 
The post-recombination annihilation signal does not depend on the (pre- 
recombination) time at which matter and antimatter domains first come into 
contact. To solve Eqs. ([T|)-(§|) we choose initial conditions at recombina- 
tion such that the matter and antimatter domains have constant density, 
have no peculiar velocity and touch along the surface x = 0. Our results 
are more conveniently presented in terms of y rather than time, accord- 
ing to dy = —yH{y)dt. For our fiducial choice of cosmological parameters 
H{y) = Hoy'/\ 

The dotted curve in Fig. ^ is the fluid temperature T{y) in a conventional 
universe. The remaining two curves are the computed temperatures of the 
B = universe: the solid curve is T{y) within the electron range where heat- 
ing by relativistic electrons dominates; the dashed curve is its value outside 
this region, where UV photons are the only heat source. For y > 400, the 
GBR is an effective heat bath keeping matter and radiation close to thermal 
equilibrium. Heating due to the annihilation products plays an important 
role at lower y: it increases the fluid temperature leading to a larger fluid 
velocity. According to Eq. (|^) the annihilation rate J is thereby enhanced. 

The solid curve in Fig. ^ is our numerical result for J, the annihilation 
rate per unit surface area defined by Eq. (^). The dashed curve is the approx- 
imation given by Eq. (H) using the temperature obtained from the numerical 
integration. Although its derivation ignored if^, Eq- ® agrees quite well 
with our numerical result for this choice of T{y). At larger redshifts, the 
motion is diffusive and our numerical result should be (and is) substantially 
less than the qualitative hydrodynamic estimate^, as is seen in the figure. 
Had we used the matter temperature of a conventional universe in Eq. (^, 
we would have obtained an annihilation rate nearly two orders of magnitude 
smaller. The heating of the fluid by annihilation debris (described by H^) has 
a dramatic effect on the annihilation rate J, and a fortiori on the consequent 
signals of the B = universe. 

At all redshifts the annihilation rate is determined by the flow of the 

^Our diffusive results for y > 400 are at variance with those of ||l^, where the annihi- 
lation rate prior to recombination is estimated on the basis of proton free-streaming. 
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Figure 2: Temperatures (in eV) as functions of redshift y = 1 + z 



highly ionized matter and antimatter fluids into the annihilation zone. The 
momentum-transfer cross section ac in proton-antiproton Coulomb colli- 
sions, which controls diffusive mixing of these fluids, is large compared to 
the annihilation cross section cXann- If mixing results only from diffusion, as 
in quasi-static laminar flow, the annihilation current would be reduced by a 
factor ~ {ca.nn/o'cY^'^ relative to J. However turbulence produces full mix- 
ing in the annihilation zone, while leaving the average flow unaffected, thus 
justifying our neglect of the Coulomb scattering term in equations (1-3). 

An analysis of fluctuations about a laminar flow into the annihilation 
zone demonstrates an instability towards turbulent mixing (this is analogous 
to the instability of a planar flame front in combustive flow [|19|.) For the 
width A of the annihilation zone we have obtained, the Reynolds number 
is TZa ~ Aacn > 10^, large enough to insure a turbulent flow at this and 
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larger scales. Turbulence efficiently mixes the fluids in the annihilation zone, 
but does not significantly retard their mean motion in the depletion zone. 
Thus turbulence drives the flow to the solution we have discussed, wherein 
the annihilation rate is determined solely by the rate at which material can 
be transported toward the annihilation zone. 

The drag on the fluid exerted by the CBR and the velocity redshift due 
to expansion suppress turbulence on large scales. The first effect dominates 
during the redshift range of interest, suppressing turbulence for scales A > 
f mp c/(crT M7), which is larger than the width A of the annihilation domain. 
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Figure 3: Annihilation rate J (in particles per cm^ per second) as a function 
of redshift y. The solid curve is our numerical solution; the dashed curve is 
an approximate result discussed in the text. 
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4 Distortion of the CBR 



Measurements of the CBR, being much more precise than those of the CDG, 
might be expected to provide the most stringent constraint on the B = 
universe. In this section, we use our conservative calculation of the annihila- 
tion rate to estimate the distortion of the CBR spectrum. In performing this 
calculation, we make several approximations that somewhat overestimate the 
effect. Nonetheless, the consequent distortion lies well below the observed 
limit, and provides no constraint at all. 

Annihilation produces relativistic electrons and energetic photons. Anni- 
hilation electrons have a direct effect on the CBR by scattering photons to 
higher energies, thereby skewing the CBR spectrum. Moreover these elec- 
trons heat the ambient plasma. The heated plasma produces an additional 
indirect spectral distortion. (The energetic photons from neutral pion de- 
cay have energies too high to have much effect on the cosmic microwave 
background.) 

To compute the direct effect, we must determine the number of CBR 
photons scattered from energy tUj to cj/ by a single electron. This function, 
(PNiujf, uji) I duf dui, is computed in Appendix |A.2| . The electron multiplicity 
per pp annihilation is similar to the photon multiplicity, measured to be 
g ~ 3.8. The number of annihilation electrons made per unit volume and 
time is gJ/d, where 1/d = y/do is the average domain surface-to- volume 
ratio at epoch y. The spectral distortion 6u^{u!) (energy per unit volume 
and energy) satisfies a transport equation: 



H{y) d{y) J \ dudoo dudoo 

We have ignored absorption of UV photons by neutral hydrogen because the 
B = universe is largely ionized. 

The direct contribution to the CBR distortion is the solution to Eq. (^ 
evaluated at the current epoch: 6u^{uj) = 6uy{uj, 1). It is given by: 

rys dy 

ivR y^ 

where we have confined the source to 1100 > y > 20, the era of unavoidable 
annihilation. To evaluate the integral we use the annihilation rate J com- 



rys dv 

5U,{U)= JLA^ujy^y)^ (10) 
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puted in Section ^. Figure ^ displays the result for a current domain size of 
20Mpc. Note that is always less than 3 x lO'^ cm'^ ~ l.gxlQ-^Tl^. 

The limit set by COBE-FIRAS on rms departures from a thermal spec- 
trum is \6uy{u) I < 7.2 X 10~^ Tq throughout the energy range Tq < a; < 10 Tq. 
This upper limit is four times larger than our computed signal for the min- 
imum domain size. Because larger domains yield proportionally smaller re- 
sults, we obtain no constraint on the B = universe. 

The indirect contribution to the CBR distortion results from a temper- 
ature difference T — between the heated ambient fluid and the CBR. It 



may be described by the Sunyaev-Zeldovich parameter Y |]22 



r= /^^^^i^^f^d/, (11) 

J rrif.C'^ 

where the integral is along the photon path dl = —cdy/yH{y). 

Within the electron range, collisions between annihilation electrons and 
the plasma result in a temperature profile T{y) shown as the solid curve in 
Fig. |[ Outside the electron range, reheating is due to photons up-scattered 
by these electrons, resulting in the temperature profile shown as the dashed 
curve. CBR photons may have traversed regions of both types. To compute 
Y, we use the higher temperature profile (the one within the electron range). 
We thereby overestimate the signal. Our result is F < 9 x 10""^, which is 
over an order of magnitude below the COBE-FIRAS limit of |y| < 
1.5 X 10^^. We conclude^ that current observations of the CBR spectrum 
yield no constraint on the B = universe. 

The energy spectrum of uplifted CBR photons shown in Fig. ^ extends 
into the visible, falling like 1/v^- Most of the energy remaining from nuclear 
annihilation resides in this tail. Nevertheless, the diffuse intensity of the night 
sky is well above this level. 



5 The Diffuse Gamma-Ray Spectrum 

In this section, we use our conservative calculation of the annihilation rate 
to determine a lower bound to the CDC signal. We find that annihilation in 
a. B = universe produces far more 7-rays than are observed. 

^An additional contribution to Y arises as CBR photons pass through transitional re- 
gions being re-ionized, but is two orders of magnitude smaller than the effect we discussed. 
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Figure 4: The CBR spectral distortion. Beyond the range shown, 6u^ oc 

1/y/uJ, up to Uj/Tq ~ 10^. 

The rehc spectrum of 7-rays consists primarily of photons from 7r° decay. 
Let ^{E) denote the inclusive photon spectrum in pp annihilation, normal- 
ized to g, the mean photon multiplicity^. The average number of photons 
made per unit volume, time and energy is J/d. These photons scat- 

ter and redshift, leading to a spectral flux of annihilation photons F{E, y) 
(number per unit time, area, energy and steradian) satisfying the transport 
equation: 

®The measured photon spectrum can be found in pC[ | and is further discussed in Ap- 
pendix 



16 



The first term on the RHS is the annihilation source and the second is a 
scattering sink. We shghtly underestimate F{E,y) by treating all scattered 
photons as effectively absorbed. In this case: 

^(^' y) = '"""Hlyf^^ ^ ' ^^^^ 

with the photon interaction cross section and ne{y) the electron density. 
For the relevant photon energies, it matters little whether photons encounter 
bound or unbound electrons. 

Integration of Eqs. (|TD-(0) gives the photon flux today, F{E) = F{E, 1): 



fVR cJ{y')^{Ey') rv ay , „ 



<-y' dy" 

If' 



(14) 



H{y')y 



/3 



Measurements of the CDG flux are shown in Fig. |^. From 2 MeV to 



10 MeV, preliminary COMPTEL satellite measurements |23| lie roughly an 
order of magnitude below|^ the earlier balloon data [^. Figure ^also shows 
our computed signal F{E). The upper curve corresponds to the smallest 
allowed domains, do = 20 Mpc, the lower curve to do = 1000 Mpc. The 
signal is linear in 1/do. The relic photon distribution is redshifted from the 
production spectrum (which peaks at ~ 70 MeV), and is slightly depleted 
at low energies by attenuation. 

Our conservative lower limit to the 7-ray signal conflicts with observations 
by several orders of magnitude and over a wide range of energies, for all 
values of do < 10^ Mpc, comparable to the size of the universe. We could 
argue that the satellite data excludes even larger domain sizes, but we would 
soon run into questions of the precise geometry and location of these nearly 
horizon-sized domains. 



6 Closing Loopholes 

Can our 'no-go theorem' for the B = universe be skirted by changing the 
input parameters, modifying our hypotheses, or including other effects? Here 
we examine the sensitivity of our conclusions to the chosen values of cosmo- 
logical parameters, to the possible existence of primordial magnetic flelds, 
and to the assumed isentropic nature of primordial density fluctuations. 

^This discrepancy is attributed by the authors of to a rigidity-dependent back- 
ground correction that the balloon experiments could not perforin. 
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Figure 5: Data |^ and expectations for the diffuse 7-ray spectrum. 



We used a flat and dark-matter-dominated universe with vanishing cos- 
mological constant. For this case, the expansion rate is given by the simple 
expression H(y) = y^^^ Hq, with Hq the Hubble constant. Other choices for 
the cosmological parameters {Qm 7^ 1 and/or Q\ 7^ 0) would alter the y 
dependence of H{y) as follows: 



dy 

y 



-H{y) dt = -Ho (1 -a)y^ + y' + 



1/2 



dt . 



(15) 



It is only through the modification of H{y) that Hq, Q^n and affect our 
results. 

We have recomputed the diffuse gamma background (CDG) for a range of 
observationally viable values of the cosmological parameters and are unable 
to suppress the signal by more than a factor of 2. The reason is easily seen. 
Equation (|12]) shows that J oc 1/H{y), and Eq. (|l|) shows that the CDG 
flux is proportional to J/H{y), and hence to H{y)~'^. To suppress the flux, 
we must increase Hijj) beyond its value at VL^ = 1, = and h = 0.75. 
No sensible value of fl/^ has much effect at y ~ 20, when most of the CDG 
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flux arises. For = 2 or h = 0.5, two borderline possibilities, the CDG 
flux would be reduced by about a factor of two, not altering our conclusions. 

We assumed that electrons produced by annihilations travel in straight 
lines. This would not be true were there primordial (or magnet ohydro dy- 
namically generated) magnetic fields in the vicinity of domain boundaries. 
Fields with sufficiently short correlation lengths and large amplitudes would 
reduce the electron range. If the magnetically-reduced range still exceeds 
D, the width of the depletion zone, the annihilation rate is increased and 
our conclusions are strengthened. If the electron range were less than D, 
electrons would deposit their energy near the annihilation zone rather than 
throughout the plasma. However, heating by UV photons alone results in 
the temperature profile plotted in Fig. |^. Because J ~ VT, the CDG signal 
cannot be reduced by more than a factor of 3 relative to our previous results. 
Thus, the existence of magnetic fields at or after recombination cannot alter 
our conclusion. 

Finally, we claimed that matter and antimatter domains must touch by 
recombination, if they are not to produce observable (and unobserved) scars 
in the GBR. Our argument depended on the absence of strictly isothermal 
fluctuations at recombination. If this hypothesis is false, matter and anti- 
matter islands could be separated by regions of vanishing baryon density, 
with a uniform photon distribution throughout. If these isothermal voids 
are so wide that they persist after recombination, annihilation might be pre- 
vented. Annihilation might also be prevented by 'wrapping' different regions 
with domain walls, whose properties are designed to block the penetration 
of thermal matter while avoiding cosmological constraints ||2^. We have not 
further pursued these contrived lines of thought. 



7 Conclusions 

Neither the notion of a universe containing islands of antimatter, nor the 
exploration of its observable consequences are new. Indeed, the literature 
includes diametrically opposed views as to the viability of such models. The 
purpose of this paper is to present a class of models (arguably, the most 
general) for which the observable universe consists of comparable numbers of 
domains containing either matter or antimatter. These models are parame- 
terized by the typical domain size today, do. Direct searches for annihilation 
radiation show that do > 20 Mpc, and future searches for antimatter among 
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cosmic rays may increase this lower bound by an order of magnitude. 

We have found constraints on a matter-antimatter universe arising from 
phenomena taking place at cosmological distances. The potentially observ- 
able signals are identified as a distortion of the CBR, and the production 
of a relic flux of diffuse gamma-rays (CDG). We have computed these sig- 
nals with conservative assumptions and considerations based on empirical 
evidence, but with as little theoretical prejudice as possible. We find that 
matter-antimatter encounters at domain boundaries are unavoidable from 
recombination to the onset of structure formation. The detailed dynamics 
underlying our calculation of the annihilation rate is complicated. The flow 
of matter into antimatter (and vice versa) is diffusive at large y and hy- 
drodynamic at low y. Furthermore, energy deposition by the annihilation 
debris plays a crucial role, increasing the annihilation rate by up to two or- 
ders of magnitude relative to what it would have been if this effect had been 
neglected. 

Part of the energy released by annihilations at cosmological distances ends 
up as microwave photons that would appear as a non-thermal correction to 
the cosmic background spectrum. However, we find that measurements of 
the CBR spectrum do not lead to a competitive constraint on the 5 = 
universe. 

High-energy photons produced by annihilations at cosmological distances 
(most of which survive to the current epoch) are redshifted to current energies 
of order 1 MeV, thereby contributing to the diffuse 7-ray spectrum. Our 
conservative estimate of the relic CDG flux far exceeds its measured value. 
Thus, we have ruled out a i? = universe with domains smaller than a size 
comparable to that of the visible universef"^ It follows that the detection of 
Z > 1 antinuclei among cosmic rays would shatter our current understanding 
of cosmology, or reveal something unforeseen in the realm of astrophysical 
objects. 
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A The Annihilation Debris 

Each pp annihilation produces g ~ 3.8 electrons and positrons, and a similar 
number of photons. The photon spectrum has been well measured ||2^ and 
may be used to infer the electron spectrum. The photon distribution peaks 
at Ey ^ 70 MeV, and the average photon energy is (E^) ~ 180 MeV. The 
mean pion energy is twice that of the photon. About 1/4 of the energy of a 
charged pion finds its way to an electron. (The muon retains ~3/4 of the 
charged-pion energy, of which ~l/3 passes to the decay electron.) Thus, we 
expect an electron spectrum peaking at Ee ~ 35 MeV with (Ee) ~ 90 MeV. 

We must determine various properties of an annihilation electron in the 
redshift interval 20 < y < 1100: its mean range, its effect on the GBR, 
the energy it deposits in matter along its trajectory, and the ionizing effect 
of its passage. Three mechanisms control the electrons' motion in the fully 
ionized plasma. With p = jS'jme and in terms of our fiducial cosmological 
parameters, they are: 

• Cosmological redshift: 



dp 



c 



R 
R 



p = Hopy'/' = Kciy)f37, 



(16) 



Kciy) 



1.3 X 10^ y^/^ 



eV 



Mpc 



• Collisions with CBR photons: 
dp 47r^ 




K,iy)P7'- 



(17) 



dt 45 

7 

K,{y) ~ 0.7 y' 



Mpc ■ 



• Collisions with ambient plasma electrons: 




dt ^ me(3'^ ylGTrriga 



) 



KM{y) 



lie J_ 



(18) 



M(y) ^ 5.5 —- . 

Mpc 
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where is the position-dependent electron number density while rioo 
is its value far enough from a domain boundary to be unaffected by 
annihilation and fluid motion. 

A.l The Range of Annihilation Electrons 

Annihilation electrons lose energy as they redshift, but this mechanism — 
given by Eq. (0) — is negligible compared with coUisional energy loss through- 
out the interval 20 < y < 1100. Collisions with CBR photons — given by 
Eq. (]T7|), for which dp/dt oc 7^ — dominate over most of the trajectory. As 
an electron becomes non-relativistic, collisions with background electrons — 
given by Eq. (]TB|), for which dp/dt oc 1//?^ — come into play. These mecha- 
nisms cross over at /3^7^ — S/y, a point denoted by Peg ileq)- Some typical 
values are peg ^ 0.62, 0.33, 0.19 at y = 20, 200, 1100. 

To compute the range ^(70, y) of an electron with initial energy •jo^e, 
we use Eq. (^) throughout its trajectory, and ignore the small effect of 
multiple-scattering corrections. For y > 20 the neglect of other energy-loss 
mechanisms leads to a negligible overestimate of L. Integrating Eq. (|l^, we 
find: ^ 

L(7o,?/) = ^ arcsin/3o ^ 0.8 X lO^M — 1 Mpc . (19) 
^7 V ^ / 

For an initially relativistic electron arcsin/?o — vr/2, and the electron range 
is insensitive to the initial electron energy. The dependence of L on 70 is 
hereafter suppressed. 

The previously established limit on domains of uniform composition is 
d{y) > 20/y Mpc. For y < 30, the electron range exceeds this minimal 
size and our one- dimensional approximation breaks down. Because we find 
a much stronger limit on the minimal domain size, this complication need 
not be faced. The result for the electron range, including all three sources 
of energy loss Eqs. (|T^-(^), is plotted in Fig. Throughout the relevant 
redshift interval, L is small compared with the horizon. 

A. 2 UV Photons 

We compute the spectral distortion caused by the passage of one electron 
(with initial energy Eq = 7o^«e) through a thermal bath of CBR photons. 
Compton scatterings conserve photon number but skew the spectrum toward 
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higher energies. The initial spectral distribution of CBR photons is dn^/duj = 
{uj/TifM^uj), with U = l/(e'^/^-^ - 1). Let d'^N{ujf,uJi)/duJidujf denote the 
number of photons transferred by one electron from the frequency interval 
dui to the interval duf. Define: 

d'^N{ujf,uJi) _ d^N{ujf,uji) dn^juj) 
duJidojf dujfdn dui 

The function d^N{ujf,uJi)/dujfdn may be regarded as the spectral distribu- 
tion of struck photons of frequency Uf produced during the voyage of one 
energetic electron through an isotropic, monochromatic photon gas of unit 
density and frequency oji. 

Let dfli{6i,(f)i) be the differential solid angle about the initial photon 
direction, and Vi be the relative speed of the colliding particles. We choose 
to measure angles relative to the total momentum direction of the colliding 
particles. The function d'^N/dufdn is obtained by averaging the differential 
transition rate over target photon directions, and integrating in time, along 
the electron trajectory: 

d'^N _ r r dfli da ^ 
dujjdn J J An ^dcuj ' 

where we have neglected the small effect of stimulated emission. 

The computation is simplified if we note that '-/T^ <^ mg, so that the 
Thompson limit applies and 

da 3aT 
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dujf 16/i4/55^io^^ 



|/i272 (1 + 2^^){1 - 27V) + (3 - 472)/iV + 

+r (r - 2 /i 7') [3 - 6 /i 7' + V (l + 2 7')] } ©( < r < ^) , 

where r = uj/uji and = 1 — /5cos6'j. 

Carrying out the integrations in Eq. (^) gives our result for d'^N/dujf dn. 
(The dt integration is most easily performed by trading dt for dp using 
Eq. (p^). This integral extends from from po — '^eTo to Peq — leg Peg- 
The result is insensitive to the y-dependence of p/.) 

A. 3 Ionization 

Here we show that the fluid is almost totally ionized by annihilation electrons 
at all relevant times. The value of the ionization fraction, x, results from a 
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compromise between the recombination and ionization rates. Annihilation 
electrons ionize the material they traverse both directly, via electron-atom 
collisions as described by Eg. (|T8|) , or indirectly, via the UV showers discussed 
in Appendix [A .21 . We discuss the latter effect, which is more important. 



Many of the photons up-scattered by annihilation electrons have ener- 
gies exceeding the hydrogen binding energy {B = 13.6 eV), and can ionize 
hydrogen atoms via 7 + if — > e + p. The photoionization cross section for 
hydrogen atoms in their ground state, ax, falls rapidly from a very large 
threshold value (Tk{B) ~ 8 x 10~^^ cm^: 

aA'(^) - <yK{B) {B/ujf Q{u - B) . (22) 

We compute the effective ionization cross section ax for the entire UV shower 
associated with a single electron by integrating the product of ax with the 
photon number distribution: 

r I 

aK = / — (^K((^f) dujf duJi ~ 1.4 X 10"^^ cm^ \/ llOO/y . (23) 

J duJi du!f ^ 

This cross section is four orders of magnitude larger than axiB) and reflects 
the large number of photons scattered by a single electron. 

The total ionization rate is the difference of the photonionization rate and 
the recombination rate. The former is obtained by multiplying the effective 
ionization cross section for a single annihilation electron ax by the flux of 
electrons. Because half of the produced in an annihilation zone move to 
either side, the flux is half the multiplicity g times the annihilation rate J. 
The total ionization rate x (per second and per baryon) is: 

X = ^ J {1 - X) Ck - rieX'^ ((Tree Ve) , (24) 

where the recombination coefficient to all states but the ground state is 
(a^ee Ve) ^ 1.14 X 10"^^ T'^^^ [l - 2.20 log T + 0.814 T^/^j cm^ s"^ . (25) 

The coefficient of 1 — x in Eq. (P^ ) is much greater than the coefficient of 
at all relevant epochs. Consequently, the ionization is very close to one: 

1 — X ~ — - — -C 1 . (26) 

gJcTK 
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In the previous argument, no allowance was made for photon absorption de- 
spite the large photoionization cross section. Because the (quasi-) equilibrium 
ionization is nearly total, UV photons are unlikely to encounter atoms. 

Near the region of electron production, the UV photon shower has not 
fully developed, so that the ionization is smaller than Eq. (^) indicates. Our 
calculation of dN/dujfdn can be modified to treat this case. We find that the 
UV flux near the annihilation zone is sufficient to maintain total ionization 
to within a few percent. 

The UV flux generated by annihilation is sufficient to prevent recombina- 
tion by producing and sustaining almost total ionization. However, for large 
values of d^, regions lying far from domain boundaries recombine as in a stan- 
dard cosmology. A moving front develops between ionized and recombined 
regions as the UV flux progresses. The velocity of the front is 

c 



where ^ is the ratio of the nucleon number density to that of the incident 
UV flux. We find Vf ~ c/3 at y = 1100 and f/ ~ c at ?/ = 20. The intense 
energy deposition taking place within the front makes an unobservably small 
contribution to the Sunyaev-Zeldovich parameter. 



A. 4 Energy Deposition 

We compute the heat function H^: the energy deposited in the plasma, per 
unit volume and time, by annihilation electrons and UV photons. 

In regions within the electron range, this function is dominated by the 
primary electron contribution. For 7 > 7eg, collisions with GBR photons 



determine the evolution of the electron velocity according to Eq. {\17\} . De- 



noting the energy deposition to matter for this portion of the trajectory by 



£1, we integrate Eq. ([T8|) to find: 

Si = KM{y) — f^' ^ = n ^27) 

Here L' is the distance traveled when j = 'jeq- 

L'iy) = {rrie/K^) arccos/3e<7 • (28) 

Most of the remaining energy, £2 — {^eq — 1) "^eC^, is deposited in matter 
over a relatively small distance interval. About one third of the energy depo- 
sition to matter takes place during this short stopping stage. In the following 
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analysis, we ignore this term, thereby underestimating electron heating by 
~ 30%, and slightly underestimating the production of CDG photons. 

Electrons arise as an isotropic flux from the thin annihilation zone of 
width A. The angular average, per electron, of the energy deposition to 
matter at a distance / ^ A from this zone is: 




(29) 



where the integration variable is the distance traveled by an electron along its 
trajectory. Within the depletion zone {dE/dl)M is a slowly- varying function 
of / that is roughly proportional to the electron density Ue- 

dE\ , eV 



/d_E^ 



(^y ^ (30) 



\ dl / M Mpc 

where 10 < a < 20. For our computations we use the smallest value of a. 

Half of the produced in an annihilation zone move to either side. Thus 
the flux is g J/2, and the electron contribution to the heat function is: 



H. = -U{-jr) ■ (31) 




The UV photon contribution is small in comparison with that of the electrons. 

Outside the electron range, only UV photons contribute to H^. In an 
ionizing collision, 'j + H ^ e + p, the mean kinetic energy 6E of the recoiling 
photoelectron is: 

1 r d^N 
Ok 



If d'^N 

6E=— uj- — — (Tie (tu)rfcjrfcu, -5-5.4 eV, (32) 
(Jk J duj duji 



where the cross sections and distribution function are those of Appendix |A.2 



The rate per unit volume of such collisions is JaK{^—x)neg/2. Using Eq. (|26D 
we express this rate as nl{a^ecVe)- Multiplying by the mean recoil energy, we 
obtain the heat function: 

= SE (cTrccfe) nl . (33) 

The UV photon flux has disappeared from this expression, reflecting the 
quasi-equilibrium state of the ionization. As a welcome consequence, is 
insensitive to additional UV photons arising from annihilation zones other 
than the nearest. 
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